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While free and weakly interacting particles are well described by a second-quantized nonlinear 
Schrodinger field, or relativistic versions of it, with various approximations, the fields of strongly 
interacting particles are governed by effective actions, whose quadratic terms are extremized by 
fractional wave equations. Their particle orbits perform universal Levy walks rather than Gaussian 
random walks with perturbations. 
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Quantum-mechanical physics is explained with high 
accuracy by Schrodinger theory. The wave equation 
for many particles can convenienty be reformulated as 
a second-quantized field theory, with an action that is 
the sum of quadratic and an interacting term 



■J\- -^int 1 

where the term A2 has typically the form 



(1) 



^2= J d^xdtij* {ji,t)[idt + h^Vy2 -V^(x)]V'(x,i), (2) 

with D being the space timension, m the mass, and V{x.) 
some external potential. The interaction term A-mt may 
be approximated in molecular systems by a fourth-order 
term in the field 



d^xd^x' dti)* {-k' , t)V'*(x, t)Fi2(x, xOV-Cx, OV-Cx', t),(3) 
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where Vi2(x, x') is some two-body potential. 

If relativistic velocities are present, the field is general- 
ized to a scalar Klein-Gordon field, or a quantized Dirac 
field. In molecular physics, the fourth-order term is due 
to the exchange of a minimally coupled quantized pho- 
ton field and is proportional to e^, where e is the electric 
charge. The field equations may be studied with any 
standard method of quantum field theory, and correc- 
tions can be derived using perturbation theory in powers 
of a = /h PS 1/137. Since a is very small, this appeoch 
is quite successful. 

If time is continued analytically to imaginary values 
t — IT, one is faced with the so-called Euclidean ver- 
sion of quantum field theory. Then perturbation the- 
ory may be understood as developing a theory of parti- 
cle physics from an expansion around Gaussian random 
walks. Indeed, the relativistic scalar free-particle propa- 
gator of mass m in Z? + 1-dimensional euclidean energy- 
momentum space p'' = (p,P4), has the form 



G{P) 



1 



ds e 



(4) 



where the energy has been continued analytically to 
P4 — —lE. The Fourier transform of e^'^'-P is the dis- 
tribution of Gaussian random walks of length s in Z? + 1 
euclidean dimensions 



P(x, X4) = {Alls) 



-(D + l)/2g-(|: 



(5) 



which makes the propagator ^ a superposition of such 
walks with lenghts distributed like e"^™'' ^M,. This 
propagator is the relativistic version of the free-field prop- 
agator of the action ([2|). The second-quantized field the- 
ory described by ([T]) accounts for grand-canical ensembles 
of orbits with their two-body interactions Q ■ 

Gaussian random walks are a natural and rather uni- 
versal starting point for many stochastic processes. For 
instance, they form the basis of the most important tool 
in the theory of financial markets, the Black-Scholes op- 
tion price theory 0] (Nobel Prize 1997), by which a 
portfolio of assets is hoped to remain steadily growing 
through hedging. In fact, the famous central-limit theo- 
rem permits us to prove that many independent random 
movements of finite variance always pile up to display a 
Gaussian distribution Q . 

However, since the last stock market crash and the still 
ongoing financial crisis it has become clear that realistic 
distributions belong to a more general universality class, 
the so-called Levy stable distribution. They are the uni- 
varsal results of a pile up of random movements of infinite 
variance 0- They account for the fact that rare events, 
which initiate crashes, are much more frequent than in 
Gaussian distributions. These are events in the so-called 
Levy tails oc l/|a;|^^'^ of the distributions, whose descrip- 
tion requires a Hamiltonian |3l 



H = const {p' 



2nA/2 



(6) 



Such tail-events are present in the self-similar distribu- 
tion of matter in the universe (sl-lioj , in velocity distribu- 



III, 



tions of many body sytems with long-range forces 
and in the distributions of windgusts |12| , oceanic moster 
waves and earthquakes 14 1, with often catastrophic 
consequences. They are a consequence of rather general 
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maximal entropy assumptions In the limit A 2, 
the Levy distributions reduce to Gaussian distributions. 

The purpose of this note is to point out, that such dis- 
tributions occur quite naturally also in many-particle sys- 
tems, provided the interactions are very strong [3l ■ They 
have been observed in numerous experiments at second- 
order phase transitions. The most accurate measurement 
of this type was done in a satellite (the so called Infrared 
Atronomical Satellite IRAS) by studying the singularity 
of the specific heat of superfluid ^He near the critical 
temperature 17[ . The observation agreed extremely well 
with the theoretical strong-coupling prediction [l^l ■ 

The field of a strongly interacting A^-body system is 
usually a multivalued function. Singularities perforate 
the space via vortex lines (for instance in type II su- 
perconductors or in superfluid ^He), or via line-like de- 
fects in the displacement field of a world-crystal formula- 
tion of Einstein(-Cartan) gravity [3l ■ If the positions of 
two particles are exchanged, one obtains a factor -1-1 for 
bosons or —1 for electrons. In two dimensions, one may 
even obtain a general phase e"^ (anyons) 22 1. 

A strongly interacting field system has a conformally 
invariant Green function |20l422| 



(7) 



If the dimension D differs only by a very small amount e 
from the critical dimension Dc, where the theory is scale- 
invariant, i.e., D = Dc+e, then 7 is of order e and z differs 
from unity by a similar amount. Such a power behavior 



is assured near Dc if the Gell-Mann-Low function [23 1 
has an infrared-stable fixed point in the renormalization 
fiow of the coupling constant. Very close to the critical 
dimension, a lowest approximation to G{p,p4) is 



G(p,P4) = {pl-'ll + DxipVpd^^']}-' 



(8) 



where A is close to 2, and Dx is a generalization of the 
diffusion constant in the Fokker-Planck equation. 

Time-independent propagators involve the limit j)4 — >■ 
0, where the correlation function behaves like 



G(p,0) oc IpI 



-2+1-1 



(9) 



The index 77 is the anomalous dimension of the field, 
which is also of order e. The existence of this limit in (|8]) 
fixes the scaling relation 



A = 2(l-7)/z = 2-77. 



(10) 



See Appendix for the calculation of the exponents to or- 
der e. The Green function ([8]) determines the probability 
distribution of particle after a time t via the double frac- 
tional Fokker-Planck equation 

[pI--' + i?A(p')^/']P(x, t) = ,5(t)5(^) (x), (11) 

where pi = dt-, p = id^ = *V. A convenient defini- 
tion of the fractional derivatives uses the same formula 



as in the dimensional continuation of Feynman diagrams 
(_V')^/2 _ r[A/2]-i/daa-V2~ieAvV2 The 

solution of ([TT|) is given in the literature [26[ and reads 



-H. 



2,1 



^D/2|x|D/2 2,3 I 2a£,^il- 



(l,l);(l-7,l-7) 



(l,l),(C/2,a/2);(l,a/2)y 



I, (12) 



where H2 '3 is a Fox H-function 27| . In the limits 7 = 
and a = 2, this reduces to the standard quantum me- 
chanical Gaussian expression {4'KDxt)~'^^^e~^^^ /'^Dxt_ 
For 7 = 0, a = 1, the result is 



P(x,0 



Dxt 



^(D+l)/2|xp+l 1-1 \ |x|2 



{l/2-D/2S)\ 



(0,1) 



(13) 



which is simply the Cauchy-Lorentz distribution function 

[T{D/2 + l/2)/7r(^+i)/2]i?;,i/[(^At)' + |xp] Wi/2. 

The probability pip may be calculated from the doubly 
fractional canonical path integral over fluctuating orbits 
t(s), x(s) ^4(5), p(s) viewed as functions of some pseudo- 
time s [28|: 

{Xbtb.sb|x„i,.s4 = j VxVtVpVp^e--^, (14) 
where A is the euclidean action of the paths t(s),x(s): 

A = j ds[i(px' - ipit') - nip, Pi)]. (15) 

Here t'{s) = dt{s)/ds, x'(s) = dx{s)/ds, and T-L{p,Pi) = 



P4 



1-7 



-\- D\{p^)'^^^ . At each s, the integrals over the com- 
ponents of p(s) and P4(s) run from —00 to 00, whereas 
those over P4{s) run from —ioo to ioo. At the end we 
obtain P(x, t) from the integral ds{xt s\0 0}. 

If 7 = 0, the path integral over P4{s) yields the func- 
tional S[t'{s) — 1], which brings (fT4)) to the canonical path 
integral 



{xbtb\Xata) = J Vx.Vpe , 



with 



J dt[ipx - Dxip 



n2^V2l 



(16) 



(17) 



Now P(x, t) = (xi|O0) satisfles the ordinary fractional 
Fokker-Planck equation 



[p4 + i?A(p')^/']P(t,x)=<5(t)<5(^)(x). 



(18) 



This has been discussed at length in recent literature [29| . 

At this place it is worth mentioning that the 
probability can be written as a superposition 
J^{da/a)fx{crt~'^^^)PQ{a,x) of Gaussian distribu- 
tions Pg(ct,x) = {AT:a)-^/'^e-'^''^'^'^ with weight 



/a(<t) = SdJ2 



(-l)"cr-"^/2 



(n-f l)!r(D-l-nA/2) 



D 



i/X 



(19) 
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where So = 27r^/^/r(D/2) is the surface of a sphere in 
D dimensions. 

If 7 ^ 0, the above functional (5-function is softened, 
and the relation between the pseudotime s and the phys- 
ical time becomes stochastic. It is governed by the prob- 
ability distribution that solves the path integral the 



{tbSb\taSa} = J 'DWp4^ CXp 1^ J ds[p4t' - pI '^)] 



.(20) 



For imaginary p4 

H{ri) which has the property that 



we define a noise Hamiltonian 



-Vi = 



drje 



-Piri-H{ri) 



(21) 



The inverse of the Fourier integral yields the noise proba- 



bility PM = 
tional 13 ll: 



, and a probability func- 



"Dpi exp 



ds {pii-j 



(22) 



Using this we may solve the stochastic differential equa- 
tion of the Langevin type 



t'{s) = i]{s), 



(23) 



in which the noise t]{s) has a zero expectation value for 
each s, and the correlation functions for n = 2,4,6,... : 

(r?(si)...r/(s2„)) = /p77 7,(si)...77(s2„)PM. (24) 



If 7 = 0, the solution of 1^ is = S[ri{s) ~ 1], 

implying that r]{s) ceases to fluctuate, and (f23| becomes 
t'{s) = 1, so that t = s. 

In the past, many nontrivial Schrodinger equations 
(for instance that of the 1 /r-potential) have been solved 
with path integral methods by re-formulating them on 
the pseudotime axis s, that is related to the time t via 
a space- dependent differential equation t'(s) = f{x{t)). 
This method, invented by Duru and Kleinert [32| to solve 
the path integral of the hydrogen atom, has recently been 
applied successfully to various Fokker-Planck equations 
33l |34| . The stochastic differential equation may be 
seen as a stochastic version of the Duru-Kleinert trans- 
formation that promises to be a useful tool to study non- 
Markovian systems. 

Certainly, the solutions of Eq. ((T8)) can also be obtained 
from a stochastic differential equation 

x = r,, (25) 

whose noise is distributed with a fractional probability 



P[ri] = / 



(26) 



Experimentally, a system with in the strong-coupling 
limit can be produced by forming a Bose-Einstein con- 
densate (BEC) in a magnetic field whose strength is 
tuned to a Feshbach resonance (ssj of the two-particle in- 
teraction. In a BEC, the four-field term in the interaction 
dS]) is local and parametrized by (x,x) cx ^^(x — x'). 
At the Feshbach resonance, the bare coupling strength 
g goes to infinity [s^, and the renormalized coupling 
f/fl, multiplied by 6/x~'^/(47r)^, converges to a fixed point 



g* K, 0.503 [see Fig. 17.1 in Ref. l21|), where /i is some 
mass scale. 

The theoretical tool to describe the physics in this 
regime is the effective action r[vl',\l'*]. This a func- 
tional of the classical expectation values of the quan- 
tum fields 4'(t,x) = ('i/)(t,x)), and contains all informa- 



tion of the full quantum theory [21|, |47| . It is the Leg- 



endre transform of the generating functional Z[rj,rj* 



of the full quantum theory, and is 



extremal on the physical field expections. All its vertex 
functions can be found from the functional derivatives 
of r[^,^'*]. In the strong-coupling limit, the effective 
interaction changes the interaction ^ to an anomalous 
power law r'"*[*,**] = {gc/2) J dtd^x |*(t, x)|'^+i. The 
power (5 is a critical exponent that is measured experi- 
mentally by the relation _B = Its value is deter- 
mined by 77 via the so-called hyperscaling relation (37l | 
6 = {D + 2 - rf) / [D - 2 + rf). The value of gc is related 
to the critical value g* w 0.503 by .9^/^"''^/^^"^+''^ = 
(25*)(''-i)/2(47r)2/24 « 6.7. As a possible application we 
may study the behavior of a triangular lattice of vortices 
which form in a rotating Bose-Einstein condensate [s^, 
and letting the magnetic field approach a Feshbach reso- 
nance. 

The results may then be compared with a calculation 
based on a new field equation that generalizes the famous 
Gross-Pitaevskii equation [s^ 



2m 



.g|*(t,x)|2 



*(t,x) = 0. 



(27) 



The new equation is obtained by extremizing the effective 
action r[1',**] = ro[^',^'*] +r'"*[*,**], where 



Fn 



dtd'^x-^\t,^)[E^-^ - Dx{f>'^f'^]^{t,i^). (28) 



By forming (5^<=*^/(5^'t(t, x), we obtain what may be 
called the fractional Gross-Pitaevskii equation: 



2\l-7)/2 '^+1 



.9c|4'(t,x)| 



*(t,x)=0. (29) 



The fractional Schrodinger equation has many prob- 
lems, such as the nonvalidity of the quantum superpo- 
sition law, the violation of unitarity of the time evolu- 
tion, and the violation of probability conservation which 
can produce nonsensical probabilities > 1 [2^. However, 
these problems exist only if we restrict ourselves to the 
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free effective action ((28)) . but this is meaningless, since 
the entire theory is only defined by the effective action in 
the strong-coupling limit — and this contains necessarily 
additional nonquadratic terms. Hence it does not possess 
free quasiparticles as in the time-honored Landau theory 
of Fermi liquids [i^ . There is always an interaction that 
invalidates the standard discussion of Schrodinger equa- 
tions. In fact, the theory of high- Tc superconductivity 
must probably be built as a true strong-coupling theory 
of this type with electrons being non- Fermi liquids (4ol | . 

The relativistic version of the entire discussion is sim- 
pler since it is based on the euclidean Green function ^ 
in which p denotes the D — 1-dimcnsional vectors (p,p4). 
The Fourier transform is the distribution fulfilling the 
Fokker-Planck equation 



and possessing the path integral representation 



(p^ 



(30) 



(31) 



The e-cxpansion is now around Dc = 4 in powers of e = 
— {D~ Dc). The critical exponent r] is small of order e^: 



0.04. It can be ignored for 6=1. 



7j = e750-H • 
power 5 in the interaction is 3 -|- e -I- 23e^/50 - 

The time-independent fractional 
equation reads now 



The 
4.76 



Gross- Pitaevskii 



(p^)^-"/^+^5ci*(x)r--i 



*(x) = 0, 



(32) 



with gc ~ 27. For a, d = D — 1 -dimensional vortex in 
D — i dimensions, it is solved by vE'(x) = a|x^|^'^ with 
A ^ (2-?7)/((5-l) = D/2-l + ri/2!ii 1/2 and for /i = 1: 
[{d + l)a^-^/4]gc = -'^CA+A-Zc^irf « 0.2, X^2~r)^. 

To compare our theory with experimental data, we 
must study the BEG in the scale-invariant strong- 
coupling limit. This is reached either by going to the 
temperature Tc of the second-order phase transition, or 
by raising the magnetic field B towards the field strength 
Be of a Feshbach resonance. Then the coherence length 
^ grows like ^ cx \t\~'^ where ly w 2/3 [U 42 1, and 
t = 1 - T/Tc or i = 1 - B/Bc ([ll). If the BEG 
is enclosed in a weak harmonic trap, this adds in the 
brackets of ((27|) a term cx |xp = R^. This is normally 
observed by the condensate density going to zero lin- 



early like 1 



= 1 - R /R^ near the border i?b (i 



the Thomas-Fermi approximation) [43[. For B near Be 
(or T near Tc), however, the anomalous power S will 
lead to the steeper approach to zero (1 — r^)^'' where 
2/3 = iy{D - 2 + ri) = 1 - 3e/10 ■ • • « 0.7, plotted in 
Fig. 1, as will be shown immediately. In addition, the 
central region is depleted. 

Let us study the appearance of a reduced mass m'^ oc 
(1 — r^) in the trap. The effective action will introduce it 




Figure 1: Condensate density from Gross-Pitaevskii equation 
(f27|) (GP,dashed) and its fractional version ((29](FGP), botii 
in Thomas-Fermi approximation where the gradients are ig- 
nored. The FGP-curve shows a marked depletion of the con- 
densate. On the right hand, a vortex is included. The zeros 
at r « Iwill be smoothened by the gradient terms in (|32p . 



into (Hal) in the form fi^-'i{m'^)''^^-'^'> f{\^\^ /{m^f^) with 
a Taylor series of f{x) (note that 1/(2 — r]) = l + ^ + -- -~ 
1.3). For small rh^, this may be resummed to a Widom 
type expression [{S + l)/Afi'!]gc\'^\^-^w(m^/\'^\^/l^) [U. 
This explains the earlier-stated steeper falloff ex. 
(?7i^)^^ of the density profiles in Fig. 1. The Widom 
function ?ii(m^/|5'|^/'^) can be expanded as 1 plus a power 
series in (rri^)"/^'' oc which contains the Wegner criti- 
cal exponent u! sa 0.8 that governs the approach to scaling 
[45| . Thereby the interaction term |^'|*~^ is modified to 
l^'I'^-^l + const X ^-"|5f|-"''/'3). Similarly, the kinetic 
terms (p^)^"**/^ in ([32|) and ([29|) approach the scaling 
limit like (p2)i-''/2[l + const xe""(p2)-"/2-f...] (4^. 

A further observable phenomenon is that the resonance 
frequency of a forced collective oscillation will depend on 
the field strength B near the Feshbach resonance [43 |. 

Summarizing we have seen that a many-body system 
with strong couplings between the constituents satisfies a 
more general form of the Schrodinger equation, in which 
the momentum and the energy appear with a power dif- 
ferent from a = 2 and 7 = 0, respectively. The asso- 
ciated Green function can be written as a path integral 
over fluctuating time and space orbits that are functions 
of some pseudotime s. This is a Markovian object, but 
non-Markovian in the physical time t that is related to 
s by a stochastic differential equation of the Langevin 
type. The particle distributions can also be obtained by 
solving a Langevin type of equation in which the noise 
has correlation functions whose probability distribution 
is specified. 

Acknowledgment: I am grateful to P. Jizba, N. Laskin, 
M. Lewenstein, A. Pelster, and M. Zwierlein for useful 
comments. 



Appendix: The lowest-order critical exponents can be 
extracted directly from the one-loop-corrected inverse 
Green function G~^(E, p) in D = 2 + e dimensions after 
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a minimal subtraction of the 1/e -pole at [48 



(33) 



For p = 0, this has a power — (— i?)^ so that 7 = ae. 
For i? = 0, on the other hand, we obtain (— p^) 
that (1 -7)/z - 1 « 7/3. 



2^1-ae/3 
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